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Forced oscillations in relativistic accretion disks and QPOs 
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Abstract. In this work we explore the idea that the high 
frequency QPOs observed in LMXBs may be explained 
as a resonant coupling between the neutron star spin and 
cpicyclic modes of accretion disk oscillations. We propose 
a new model for these QPOs based on forced oscillations 
induced in the accretion disk due to a stellar asymmetric 
rotating gravitational or magnetic field. It is shown that 
particles evolving in a rotating non-axisymmetric field are 
subject to three kinds of resonances: a corotation reso- 
nance, a Lindblad resonance due to a driving force, and 
a parametric resonance due to the time varying cpicyclic 
frequencies. These results are extends by means of 2D nu- 
merical simulations of a simplified version of the accretion 
disk. The simulations are performed for the Newtonian 
gravitational potential, as well as for a pseudo-general rel- 
ativistic potential, which enables us to explore the behav- 
ior of the resonances around both rotating neutron stars 
and black holes. Density perturbations are only significant 
in the region located close to the inner edge of the disk 
near the ISCO where the gravitational or magnetic per- 
turbation is maximal. It is argued that the nearly periodic 
motion induced in the disk will produce high quality factor 
QPOs. 

Finally, applying this model to a typical neutron star, 
we found that the strongest response occurs when the fre- 
quency difference of the two modes equals either the spin 
frequency (for "slow rotators" ) or half of it (for "fast ro- 
tators" ) . The two main excited modes may both be con- 
nected to vertical oscillations of the disk. We emphasize 
that strong gravity is not needed to excite the modes. 



1. Introduction 

Quasi-periodic oscillations (QPOs) have been observed in 
accretion disks around neutron stars, black holes, and 
white dwarf binaries with frequencies ranging from a 
few 0.1 Hz up to 1300 Hz. Recent observations have shown 
a strong correlation between the low and high frequency 
QPOs (Mauche [20021 Psaltis et al. IT999)) . This relation 
holds over more than 6 orders of magnitude in frequency 
and strongly supports the idea that the QPO phenomenon 



is a universal physical process independent of the nature 
of the compact object. 

To date, quasi-periodic oscillations (QPOs) have been 
observed in about twenty Low Mass X-ray Binaries 
(LMXBs) containing an accreting neutron star. Among 
these systems, the high-frequency QPOs (kHz-QPOs) 
which mainly show up in pairs, denoted by frequencies 
v\ and 1/2 > v\, possess strong similarities in their fre- 
quencies, ranging from 300 Hz to about 1300 Hz, as well 



as in their shapes (van der Klis 12000]) . 

Several models have been proposed to explain the kHz- 
QPOs in LMXBs. A beat-frequency model was introduced 
to explain the commcnsurability between the twin kHz- 
QPOs frequency difference and the neutron star rota- 
tion. This interaction between the orbital motion and the 
star rotation happens at some preferred radius. Alpar & 
Shaham (1985]) and Shaham ()1987p proposed the mag- 
netospheric radius to be the preferred radius. The sonic- 
point beat-frequency model was suggested by Miller et 
al. (|1998p . In this model, the preferred radius is the point 
where the radial inflow becomes supersonic. 

The relativistic precession model introduced by Stella 
& Vietri (|1998| I1999|) makes use of the motion of a sin- 
gle particle in the Kerr-spacctimc. In this model, the 
kHz-QPOs frequency difference is related to the relativis- 
tic periastron precession of weakly elliptic orbits while 
the low-frequencies QPOs arc interpreted as a conse- 
quence of the Lense-Thirring precession. Abramowicz & 
Kluzniak (|200ip introduced a resonance between orbital 
and epicyclic motion that can account for the 3:2 ratio 
around Kerr black holes leading to an estimate of their 
mass and spin. The 3:2 ratio of black hole QPOs frequen- 
cies is well established (McClintock & Remillard, [2003). 
In other models, the QPOs are identified with gravity or 
pressure oscillation modes in the accretion disk (Titarchuk 
et al. IT9981 Wagoner et al. l200i"j) . Rezzolla et al. (f2003|) 
suggested that the high frequency QPOs in black hole bi- 
naries are related to p-mode oscillations in a non Keplerian 
torus. 

We propose a new explanation of this phenomenon 
based on a resonance in the star-disk system arising 
from the response of the accretion disk to cither a non- 
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axisymmctric rotating gravitational field, for a hydrody- 
namical disk, or to a non-axisymmetric rotating magnetic 
field, for an MHD disk. 



2. Hydrodynamical disk 

We first consider a hydrodynamical disk evolving in an 
asymmetric gravitation field imposed by the rotating ac- 
creting compact object. A simplified linear analysis reveals 
the main characteristic of the perturbed motion in the 
disk. This is then confirmed by full non-linear numerical 
simulations. 



2.1. Linear analysis 

By perturbing the set of hydrodynamical equations gov- 
erning the evolution in the disk around the equatorial 
plane, and introducing the Lagrangian displacement £, 
the weak oscillations in the radial and vertical directions 
can be cast into a partial differential equation for the ra- 
dial and vertical displacement respectively. Neglecting the 
sound propagation term (not important because not lead- 
ing to resonance conditions), the Laplacian displacement 
£ r / z satisfies : 



£> 2 6 9 , id,,.. 

~W + + Tr ~dr (rP?r) 9l = 91 



p dz 



(1) 
(2) 



We have introduced the radial/ vertical epicyclic fre- 
quency k t / z and the perturbation in radial/ vertical gravi- 
tational field by Sg r / Z . The convective derivative is denoted 
by D/Dt = dt + Qdp where ft is the unperturbed local or- 
bital frequency in the disk. The 3 terms in Eq. |T]) and (J2|) 
are: a harmonic oscillator at frequency K r / Z , second term 
on the left hand side, with a periodically perturbed eigen- 
frequency, third term on the LHS, and a periodic driven 
source, on the right hand side. 

A careful analysis of this equation shows the emergence 
of three kinds of resonance corresponding to : 



2k t / z /h, 



— a corotation resonance ; 

— an inner and outer Lindblad resonance ; 

— a parametric resonance when m |f2* — fi| 
where f2* is the stellar spin, K r / Z the radial and vertical 
epicyclic frequencies, n is a natural integer and m the 
azimuthal mode of the gravitational perturbation. 

For a Newtonian disk, we have Q = n r = k z , and the 
parametric resonance condition simplifies into : 



m±2/n 



(3) 



As a consequence, the resonances are all located in the 
frequency range SI € [fi*/3, 3 Q*]. 



For the general-relativistic disk, it splits into the two 
following cases : 



2/c r/z (r,q) 
s2(r, a) ± ■ = 



m n 



(4) 



For a typical neutron star, we choose: 

— mass M* = 1.4 M ; 

— angular velocity — Sl*/27r — 300 — 600 Hz ; 

— moment of inertia /* = 10 38 kgm 2 ; 

— angular momentum a* = c/„ £l*/G M^:. 

The angular momentum is then given by a» = 5.79 * 
10~ 5 51, . For the selected spin rate of the star we find a* = 
0.109 — 0.218 and so the vertical epicyclic frequency is close 
to the orbital one k z « £1. Thus for the vertical resonance, 
we are still close to the Newtonian case given by Eq.Q. 

The results for an accretion disk evolving in a New- 
tonian and a Kerr spacetime are shown in Table [1] and [2] 
(Petri [2006]). 





Orbital frequency v(r, a) (Hz) 


Mode 


v« = 600 Hz 


v« = 300 Hz 


m 


n = 1 n = 2 


n — 1 n = 2 


1 


-600 / 200 — / 300 


-300 / 100 - / 150 


2 


— / 300 1200 / 400 


- / 150 600 / 200 


3 


1800 / 360 900 / 450 


900 / 180 450 / 225 



Table 1. Orbital frequencies at the parametric vertical 
resonance for the first three order n in the case of a New- 
tonian gravitational potential. The value on the left of the 
symbol / corresponds to the absolute value sign taken to 
be - and on the right to be +. 



2.2. Two-dimensional simulations 

From the analytical analysis of the linear response of a 
thin accretion disk in the 2D limit, we know that waves are 
launched at the aforementioned resonance loci. They prop- 
agate in some permitted regions inside the disk, according 
to the dispersion relation obtained by a WKB analysis. We 
confirm and extend these results by performing non lin- 
ear hydrodynamical numerical simulations using a pseudo- 
spectral code solving Euler's equations in a 2D cylindri- 
cal coordinate frame. Simulations were performed for the 





Orbital frequency z/(r, a*) (Hz) 


Mode 


v, = 600 Hz 


v* = 300 Hz 


771 


n — 1 n — 2 


n — 1 n — 2 


1 


— / 200 — / 300 


- / 100 - / 150 


2 


— / 300 1198 / 400 


- / 150 599 / 200 


3 


1790 / 360 899 / 450 


898 / 180 450 / 225 



Table 2. Same as Tab. [T] but for the general-relativistic 
disk. 
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t=68 . 9316T, 
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Fig. 1. Final snapshot of the density perturbation Sp/po 
in the accretion disk evolving in a quadrupolar perturbed 
Newtonian potential. Time is normalized to the spin pe- 
riod. The m = 2 structure emerges in relation with 
the m = 2 quadrupolar potential perturbation. 



t/T.=61.3491 




-60 -40 -20 20 40 60 



Fig. 2. Final snapshot of the density perturbation in the 
accretion disk evolving in a perturbed Newtonian poten- 
tial having many azimuthal modes m. 



Newtonian as well as for a pseudo-Newtonian potential. 
For instance, in the Newtonian potential, for a perturba- 
tion of mode m — 2, the stationary density perturbation 
is shown in Fig. [TJ The forbidden region where no wave 
propagation is possible, is located between r = 24.0 and 
r = 49.0 as predicted by the WKB analysis. Outside this 
region, an m = 2 spiral structure forms and rotates at the 
star speed. The corotation resonance located at r = 40.0 
is not seen at this stage. When many azimuthal modes are 
excited at the same time as in Fig. [2J it becomes hard to 
distinguish between the different permitted and forbidden 
regions for each m. No clear symmetric pattern emerges 
from the simulations due to overlapping of many m. Nev- 
ertheless, in all case, Newtonian or pscudo-Schwarzschild 
and pseudo-Kerr, the simulations agree well with the de- 
tailed linear analysis. 



vertical directions, the Laplacian displacement £ r / z satis- 



fies : 




D% 
Dt 2 


K * dz 


D% 
Dt 2 


- 2 d 

K ' r d'r 



MOP 

MoP 



dB r 
dz 

Or 



d f B r SB 
dz \ 
_d_ 
Or 



Mop 
B Z SB 



(5) 



MoP 



Now, the perturbation in magnetic field 8B r / z replaces 
the perturbation in gravitational field. We recognise again 
a harmonic oscillator with periodically varying eigenfrc- 
quency superposed on a driven source term. Therefore, 
despite the presence of a magnetic field, the behavior of 
the accretion disk, at least in the linear stage of its evo- 
lution, agrees with the study already made for a hydro- 
dynamical disk. The orbital frequencies for resonance are 
therefore also given by tables [1] and [5] for a Newtonian and 
a relativistic disk respectively. 



3. MHD disk 

We extend the previous idea to a magnetised accretion 
disk evolving in the rotating asymmetric magnetic field 
field of the accreting compact object (linear analysis and 
numerical simulations) (Petri I2005ap . This case is well 
suited for binaries containing a neutron star or even for 
white dwarfs in cataclysmic variables for which the QPO- 
phcnomcnology seems to be identical (Warner et al. l2003|l . 

3.1. Linear analysis 

Here again, applying the linear perturbation theory to the 
set of MHD equations governing the motion in the disk 
and focusing only on resonance terms in the radial and 



3.2. Two-dimensional simulations 

We again performed 2D numerical simulations by solv- 
ing the magnctohydrodynamical equations for the accre- 
tion disk. This is done by extending the previous pseudo- 
spectral method to a simplified version of the 2D MHD 
accretion disk. 

Before the time t = 0, the disk stays in its axisymmet- 
ric equilibrium state and possesses only azimuthal motion. 
At t = 0, we switch on the perturbation by adding an 
asymmetric rotating component to the magnetic field. We 
then let the system evolve during more than one thousand 
orbital revolutions of the inner edge of the disk. 

We ran a simulation in which the rotation of the star 
is taken into account. This shifts the location of the ISCO 
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Fig. 3. Final snapshot of the density perturbation in the 
accretion disk evolving in a perturbed pseudo-Kerr poten- 
tial with a = 0.5. The outer Lindblad resonance is not on 
the grid. 



closer to the surface of the neutron star as compared to 
the non-rotating case. 

We chose a star with an angular momentum of a* = 
0.5. Therefore, the disk inner boundary corresponding 
to the marginally stable circular orbit approaches the 
horizon. An example of the density perturbation in the 
pseudo-Kerr metric is shown in Fig.[3]for the m = 2 mode. 
The inner Lindblad radius is clearly identified while the 
outer Lindblad radius is outside the simulation box. 

4. Slow vs fast rotator 

Focusing on accreting neutron stars in LMXBs, observa- 
tions reveal that they can be divided into two categories 
(van der Klis I2004|) : the slow rotators possessing a ro- 
tation rate z/* w 300 Hz, and for which the frequency 
difference between the two peaks is around Av ps and 
V2 ~ 2>V\j2 and the fast rotators having i/» « 600 Hz, for 
which this difference is around Av « v*j1. 

The model presented in this work can account for this 
segregation if the innermost stable circular orbit (ISCO) is 
taken into account. Indeed, for a typical neutron star, the 
orbital frequency at the ISCO is v\ sco = 1571 Hz which is 
therefore the upper limit for any QPO frequency Discard- 
ing the resonance frequencies in the relativistic disk which 
are higher than z^ sco , we conclude from Table [2] that : 

— for slow rotators, the two highest frequencies are less 
than fj sco and given by v\ = 599 Hz and v 2 = 898 Hz, 
therefore Av = 299 Hz which is very close to . More- 
over vi w 3 in accordance with observations ; 



— for fast rotators, the highest frequency is not observed 
because the resonance is located inside the radius 
of the ISCO (zyi sco < 1790 Hz). Therefore the two 
highest observable frequencies are v\ = 899 Hz and 
V2 = 1198 Hz, having a difference Av = 299 Hz which 
is close to v*/2. 

As already claimed in the previous section, these conclu- 
sions apply to hydrodynamical (Petri I2005b|) as well as 
to magnetised (Petri [2"005cp accretion disks. Observations 
therefore strongly support our resonance model. 

5. Conclusion 

The consequences of a weakly rotating asymmetric stel- 
lar gravitational or magnetic field on the evolution of a 
thin accretion disk are as follows. Corotation, driven and 
parametric resonances are excited at some preferred radii. 
The kHz-QPOs are interpreted as the orbital frequency 
of the disk at locations where the vertical response to 
the resonances is maximal. The 3:2 ratio is predicted for 
the strongest modes and a clear distinction exists between 
slow and fast rotators, a direct consequence of the pres- 
ence of an ISCO. Nevertheless general relativistic effects 
are not required to excite these resonances. They behave 
identically in the Newtonian as well as in the Kerr field. 
Therefore the QPO phenomenology is explained by the 
same picture, irrespective of the nature of the compact 
object (black hole, neutron star or white dwarf). Indeed 
the presence or the absence of a solid surface, a magnetic 
field or an event horizon plays no relevant role in the pro- 
duction of the X-ray variability. 
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